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Abstract 

Universe evolution, as described by Friedmann's equations, is deter- 
mined by source terms fixed by the choice of pressure x energy-density 
equations of state p (p). The usual approach in Cosmology considers equa- 
tions of state accounting only for kinematic terms, ignoring the contribu- 
tion from the interactions between the particles constituting the source 
fluid. In this work the importance of these neglected terms is emphasized. 
A systematic method, based on the Statistical Mechanics of real fluids, 
is proposed to include them. A toy-model is presented which shows how 
such interaction terms can engender signiflcant cosmological effects. 
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1 Introduction 

The greatest achievement of Physical Cosmology is, up to now, the Big Bang 
standard model. This model is based on two principles which are consistent with 
large-scale observations: the cosmological principle and the universal time prin- 
ciple. The first postulates that the universe space-section F is homogeneous and 
isotropic; the latter states that the topology of the space-time four-dimensional 
manifold M is the direct product A4 = M x F, and results from the adoption of a 
cosmological time as parameter of the manifold foliation [TJ [5] . These principles 
lead [31 [3] to the Fricdniann-Lemaitrc- Robertson- Walker (FLRW) intervajj 
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^ To whom all correspondence should be addressed. E-mail: leo@ift.unesp.br. 
^ We shall be using natural units, c = 1, ?i = 1. 
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On the other hand, the source terms are described by the perfect fluid energy- 
momentum tensor. This tensor and the above interval, when substituted 
into the Einstein's equations of General Relativity 

n 1 A SttG 

(with a cosmological-constant A-tcrm) lead to the Fricdmann equations for the 
scale factor a [t) : 

b? A SttG fc 

a A 47rG , „ , 

- = ^-^(P + 3p) . 4 
a 6 6 

The solution for a (t) is obtained after inserting in these equations the consti- 
tutive relation between the pressure p and the energy-density p of the cosmic 
source fluid. 

The relation between p and p is given directly by an equation of state (EOS) 
of type p = p{p), or indirectly via the distribution functions. Either way, here 
are implicitly in use the methods of Statistical Mechanics (SM), which can, 
for pedagogical purposes, be divided into two branches [5]: equilibrium SM, 
with its ensemble technique and its partition functions; and non-equilibrium 
SM, whose distribution function entail, for example, kinetic equations as the 
Boltzmann and the Vlassov-Landau equations. The distribution function are 
not only useful to describe systems out of equilibrium, but also systems in 
thermodynamical equilibrium, to which a temperature can be attributed. The 
ensemble approach, however, can only be applied to systems in equilibrium. 

Cosmology uses both sections of Statistical Mechanics, e.g.: the EOS ob- 
tained from a partition function leads to the description of the radiation- 
dominated era in the thermal history of the universe and leads also to the 
ACDM model for the present-day cosmic dynamics [SI \7\ ; the distribution func- 
tion determined from Boltzmann equations are used in the perturbative cos- 
mological models describing the cosmic microwave background anisotropies, as 
well as the formation of structures [5]. But the fact is that, perturbative or 
not, the present formulation of Cosmology does not take into account 
the interactions between the constituents as direct sources of gravi- 
tation. Only kinematical terms are computed in the ordinary (non-perturbed) 
cosmological models, interactions being just used to explain pair-production and 
thermalization; the dynamical terms -which could be introduced via particlc- 
to-particlc potentials or through the S-matrix - arc simply forgotten. Even the 
perturbative approach [8] considers the interaction so insufficiently that the 
pressure function appearing in the expression of the energy-momentum tensor 
presents kinematical terms solely (3]. 

In section [21 we shall make explicit the absence of dynamical terms (interac- 
tion terms) both in the conventional cosmological models developed under the 
hypothesis of thermodynamical equilibrium - with their typical EOS - and in 
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the perturbative models (out of equilibrium) ~ built on the distribution func- 
tions. Based on the Statistical Mechanics of systems with interactions [5l [10] 
we discuss, in section [31 how to include those terms in Cosmology. Section [Jj 
presents a toy-model which shows the decisive influence that interaction terms 
can have on cosmic evolution. The conclusions arc summed up in section [S] 

2 Absence of interaction terms 

This section deals, in the context of the standard model, with an ideal source 
cosmic fluid, interactions between components being supposed absent. We leave 
to the next section the discussion about the real need for the inclusion of inter- 
action terms in some phases of the universe evolution. 

First, it is important to establish the period during which equilibrium Sta- 
tistical Mechanics is applicable. This is restricted by the fact that the universe, 
described by the FRW metric ([!]), is growing up — the volume containing the 
cosmic fluid is expanding. In the primeval universe (kT > 20 MeV, where k is 
the Boltzmann constant), the typical reaction rates Tpri involving the different 
constituents are much larger than the expansion rate Hpri i.e. 

^ pri ^ ^pri — ■ y'J ) 

For the energy values then prevalent the thermodynamic notion of quasi-static 
expansion holds: for each infinitesimal variation of volume, all constituents 
of the fluid are at the same temperature — they keep themselves in thermal 
equilibrium. Notice in advance that, as the curvature is neglig ible [k = Ojl 
in the early universe, the ordinary Statistical Mechanics defined on euclidean 
3-space E'^ can be used. We emphasize that Eq.fO is valid even during an 
accelerated expansion. 

The components of the cosmic fiuid will decouple progressively as the tem- 
perature decreases [111 15]. A natural question turns up: once decoupled from 
the other fluid components, will a given particular component remain in ther- 
mal equilibrium ? This question is better formulated in terms of the component 
distribution function: is there an equilibrium distribution function for the de- 
coupled component valid in general in an expanding plane space ? The answer 
is no [12]. The proof of this statement is not trivial, but it is related to the fact 
that we do not have spatially constant time-like Killing vectors in plane FRW 
space-times [9]. 

Nevertheless, it is possible to construct distribution functions for the ex- 
panding space in particular cases, as in the presence of non-relativistic and 
ultrarclativistic components. This is comforting enough since, from the early 
to the present-day universe, the decoupling stable particles are either ultrarcl- 
ativistic - photons and neutrinos - or non-rclativistic (baryonic and dark mat- 
ter)|l We are therefore allowed to use equilibrium SM when working with this 

In agreement with the recent observational data, we will assume fc = all along the text. 
^ We adopt cold (non-rclativistic) dark matter from the bcgining. 
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content. The non-equilibrium treatment will be required only for perturbative 
approximations . 

Next sub-section deals specifically with the description given by the cosmo- 
logical models in thermodynamical equilibrium, and 12. 21 discuss the models with 
components out of equilibrium. 



2.1 Cosmology for systems in equilibrium 

The classical Cosmology text-books, e.g. [U [HI [HI E] , teach us how to determine 
the pressure pi, the energy density pi and the numerical density n.i of the i-th 
component of the fluid in thermal equilibrium: 

f fMd^p, (6a) 



(27r) 

/ E,{p)Mp)d^P, (6b) 



(2^) 

9^ f IpI^ 



" (2ir)' J -SE.ip)- 

where gi is the degeneracy degree; Ei, the dispersion relation 

Ef = + mj ; (7) 
and fi{p) the distribution functions given by 

The lower (upper) sign refers to the Fermi-Dirac (Bosc-Einstein) statistics. 

The sum over all the i components is taken into the Friedmann equations 
([21 [Ij). Standard texts on Statistical Mechanics [TD] derive Pi, pi and in 
the ensemble formalism for ideal relativistic quantum particles. The grand- 
canonical partition function or the potential fJ^ for the i-th component are 
given as 

17,(F,/3,M.) ^ i^i^^^Mli^ = ± I In [1 ± .,e-n d^p = 



(27r) 

(27r)3 A. 



h^^~'''''d'P^ (9) 



where 



= e^'^' = eKA^f'+'-O (10) 

is the fugacity. In ([TOj), the chemical potential = fj,^^ + rrii has a non- 
relativistic contribution fj.^^, which is the usual term appearing in text-books. 
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We have above included also the term rrii associated to the rest-energy of the 
particle under consideration. The quantities pi, pi and rii are, then, 

Pi = ; (11a) 

P 



■ (lib) 

V, p 

(11c) 



da, 

dzi 

These prescriptions arc equivalent to the definitions ([6aj) . (j6b[) and ([6c|) . 

Analyzing the ultra-relativistic limit, kT ^ m, in Eas. (|6b[l6cp and assuming 
kT ^ fi (as is the case for the photons), we obtain 

P.-^, (12) 

the familiar equation of state for radiation. 

Conversely, if wc take the non-rclativistic limit kT <C m in the same Eas. (j6bi 
[Sc| and neglect quantum effects, it results 

Pnr = riNRm , (13a) 
Pnr = nMRkT . (13b) 

The quantum effects will be negligible if the condition nX^ ^ 1 is satisfied, A 
being the thermal wavelength 



This will hold when no particle invades any other's effective volume, of which a 
rough estimate is A'^: this is the meaning of the condition above. 

Pressure p only appears in the Fricdmann equations added to p. If we 
consider kT <C m, Eqs. (|13al Il3bp will say that 

PM = (15) 

is a very good approximation for the equation of state for non-rclativistic matter 
M (the dust approximation). 

The Fricdmann equations, together with (fT2|) and (fT5|) . enable us to calcu- 
late simplified, but analytic, solutions for the universe evolution determined by 
a{t) [3- Present-day universe, for example, is well described by the solution ob- 
tained after inserting EOS in ([21 H]), keeping A non-null, and neglecting the 
contribution of radiation — the so-called ACDM model. For the early universe, 
period during which radiation dominates, the suitable EOS is ([T^ . 

Anyway, neither the EOS simpler forms (|12[ [T5)) nor the complete expres- 
sions, Eas. (|llal lllbl lllcp . take interactions between the constituents into ac- 
count. Partition function ([9]) includes all the possible states of the z-th com- 
ponent in an ideal gas, of non-interacting fermionic (or bosonic) relativistic 
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particles. The first term of series ^ refers to the classical description (Boltz- 
mann statistics), while the others are order-by-order quantum corrections. We 
shall see in section |3] that the expression of the partition function in terms of a 
series suggests a mechanism of inclusion of the interactions order by order. 

We remind the reader that our interest here is to analyze EOS based on first 
principles; Eqs. p^ and (|13bp are examples of this class. We shall exclude any 
EOS used in the context of Cosmology constructed on phcnomcnological basis, 
as EOS for scalar fields [H [H [H [H [17] , EOS for the Chaplygin gas [Hill 
HOI mi 112 and Van der Walls' EOS [13], to mention but a few. 



2.2 Cosmology for systems out of equilibrium 

As said before, the use of non-equilibrium Statistical Mechanics is essentially 
required only in pcrturbative cosmology. According to [U [9] , perturbations on 
the FRW metric come from interactions between the components of the cosmic 
fluid as described by the Boltzmann equation, 

f=T.Cm, (16) 

j 

where fi is the distribution function of the i-th component and Cij[fi] is the 
collision term for the («, j) pair of components. It is written as 

a,[/,:(pi)]= mP3)fM) - Mpi)W2)] , (17) 

P2, P3, Pi 

where Mij is the scattering amplitude of two interacting particles with incoming 
momenta pi and P2 and out-coming momenta and ^4. 

Once the distribution function is evaluated (via Boltzmann equation) for all 
components of the fluid, the energy-momentum tensor is obtained as 

E ».//.(.-) ^iS.. (-) 

where p^^ is the comoving momentum encapsulating, of course, the energy and 
the tri-momentum p, 

Po=E; p2 ^ 5„pV . (19) 

As sources in Einstein's equations, T"" is the energy density ()6bp and T" is 
the pressure ((6c|). At first sight one could think that the Boltzmann equation 
(|16p introduces interaction terms as direct sources of gravitation, but in fact 
the usual pressure expression only partially considers dynamical (interaction) 
effects. Interactions are taken into account only indirectly, through deformations 
in the distribution functions fi and fj. Nevertheless, the distributions used are 
those of free particles. An example of dynamical pressure appears in the van 
der Walls equation 

nkT „ , , 

^-(1-4^-^'^' 
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in which A and B [3] are constants phcnomcnologically fitted for each gas. It 
cannot be obtained from ([18]). 

Let us see in more detail what happens. The deduction of the Boltzmann 
equation uses the fact that the number of pairs of particles with velocities Vi 
and Vj at time t is 

Mxi,v^)fj{xj,Vj) . (21) 

On the other hand, Statistical Mechanics tells us that this is true only when 
the correlations between particles are negligible. In fact, the number of pairs of 
particles with velocities Vi and Vj is determined by the two-particle distribution 
function 

f.ij{xi,Xj,Vi,Vj) (22) 

and, in general, 

fij{xi,Xj,v^,Vj) ^ fi{x.i,Vi)fj{xj,Vj) . (23) 
The A^-particle distribution function /Ar(xi, xat) may be expressed as [5] 

N 

fN{xi,...,XN) = Y\_M^^) + 9n{xi, ■■■,xn) ; (24) 

i 

qn determines the complete correlation degree of the system. For practical 
reasons, it is convenient to split the set of N particles in all the possible disjoint 
subsets containing at least one particle, i.e. 

fi{Xi) = fi{Xi) , 

fij {xi, Xj) fi {xi ) fj {xj ) + gij (xj , Xj ) , 

fijk{xt,Xj,Xk) = .fi{xi)fj{xj)fk{xk) + ft{xi)gjk{xj,Xk) + (25) 
+fj{xj)gik{xt,Xk) + fk{xk)9ij{xi,Xj) + gijk{xt,Xj,Xk) , 

etc. In these equations /; is the one-particle distribution function; fij is the 
two-particle distribution function, and so on. 

This is the cluster expansion formalism, function gs{xi, ...^Xs) being the ir- 
reductible correlation function of s particles. These are the functions describing 
interactions. For example, the van der Walls equation can be obtained from the 
two-particle correlation gij{xi,Xj) and a couple of suitable approximations. 

Notice that, in spite of fi being always positive, the function fi-^^...iN is not 
necessarily so. This fact is exemplified by Ea. ((20|) . where the relative values of 
A and B determines if the pressure is positive or negative!! 

A remark: the attentive reader will have noticed that in the passage from 
([23]) to ([24]) the dependence on Vi has been supressed. That is because the 
correlation functions usually depends only on the position. The /at depend both 
on position and velocity, but the sector correspondent to the velocities may be 
separated, resulting in non-correlated distribution functions fi(vi). Figure 1 
shows schematically the correspondence between the distribution functions and 
the clusters for a three-particle system. 



We are not considering phase-transitions effects. 
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Figure 1 : Cluster schematic representation for a three-particle system. 

The distribution functions for a A^-particle system are determined, order by 
order, by a hierarchical set of equations known as BBGKY - after Bogoliubov, 
Born, Green, Kirkwood and Yvon [24j. The Boltzmann equation is itself an ap- 
proximation of the BBGKY system reducing this set of linear coupled equations 
to just one non-linear equation. Considering interaction in this equation will 
solely modify the one-particle distribution functions fi. The statistical systems 
modeled by Boltzmann equation do not exhibit correlation terms in the dis- 
tribution functions. Therefore, the thcrmodynamical quantities, like pressure, 
do not present dynamical terms. This leads to the conclusion that the usual 
approach of pcrturbative cosmology does not compute interaction processes as 
direct sources of gravitation. 

We shall below discuss in which periods of cosmic history interaction should 
effectively contribute to gravitation, and how to include it. 

3 On the inclusion of interaction terms 

Qualitatively, the cosmic fluid can be separated into three different sectors: a 
kinetic part T, associated to velocities; a dynamical component V, related to 
interactions; and a massive contribution M coming from the rest mass of the 
constituents. The condition for the interaction being a relevant source is 

V{At) > T{At) + M{At), (26) 

where At indicates a given cosmological period. 

As far as interactions are concerned, each period presents its own charac- 
teristics. We will deal with two of then: the pre-nucleosynthesis cosmological 
universe (PNU) - kT > 20 MeV or red-shift z > 10^^ - and the recent universe 
- red-shift z < 20. 

In the PNU, the high temperature of the fluid warrants the existence of a 
large variety of particles, such as 7, v, tt, most of then interacting mu tuallyll 
The answer to the question about the relevance of interactions for the primeval 
cosmology comes from an involved analysis of the interacting fluid as a whole: 

^ The magnitude and type of interaction depends, of course, on the espocies considered. 
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the interaction terms may be positive or negative, depending on the nature of 
the interaction, and global cancellation may possibly turn up. 

Recent data from RHIC (Relativistic Heavy Ion Collider) |25| indicates that 
a fluid at very high temperatures (few hundreds of MeV) presents a strong 
interaction between its constituents, even possibly generating a "liquid" state 
for the hadronic matter, the CGC {Color Class Condensate) [26| . Quantum 
Chromodynamics suggests that this system is constituted basically by three 
quarks (u, d, s) and gluons. In this energy range - which corresponds to the PNU 
- the strong interaction between quarks (generically, q) overcome their kinetic 
and rest energies. We have, so, Vq{AtpNu) > Tq{I^tpNu) + Mq{I^tpNu)- Even 
then, the primordial cosmic fluid has others species that should be considered 
when applying criterion (|26p . In section |4] we exhibit a very simplified model, 
introducing interaction terms, to describe the pre-nucleosynthesis cosmological 
universe. 

Among the various components of the recent universe, one of the most im- 
portant for its evolution is non-relativistic matter. Indeed, at the present-day 
time to (red-shift z = 0) , the rest mass of this component Matjj corresponds to 
about 30% of the universes' total content [27|, and, as we go back in time, it 
becomes more and more relevant j28j . Non-relativistic matter responds to the 
Newtonian gravitational interaction Vf^ji responsible for structure formation 
and evolution, and the experimental data indicates the influence of the poten- 
tial is more effective as z diminishes (i.e., the structures grow) [52]. As the 
kinetic term T/v/j is negligible compared to M^p, the importance of the gravi- 
tational interaction Vnr is measured by its direct comparison with Mmr- The 
Newtonian potential is a long range interaction and one could ask whether it 
can decisively contribute with dynamical terms that could influence the present 
day cosmic behavior. Wc will not study this subject in the present work, but 
only mention three effects controlling its importance: 

1. As the scale factor a{t) grows, so does the interaction distance d = a{t)r 
(where r is the comoving distance) , reducing the cosmological contribution 
of Vnr- 

2. During a decelerated (accelerated) expansion, the comoving horizon in- 
creases (decreases) consequently increasing (decreasing) the global effects 
of Vnr- 

3. The basic constituents of the universe are different at each phase of the 
universe evolution, starting with fundamental particles (nuclcons, etc.), 
passing to hydrogen clouds and then to galaxies and clusters. 

Let us return to the analysis of the pre-nucleosynthesis cosmological period. 
Next section presents a prescription to include interaction in that period. 

3.1 Equation of state with interaction terms 

The inclusion of interactions via EOS can be done by the ensemble formalism 
through the pcrturbative treatment of real gases. Mayer and collaborators devel- 
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oped in 1937 the systematic approach of cluster expansion for a non-relativistic 
classical (non-quantum) system [30]. Just after that, in 1938, Kahn and Uh- 
lenbeck began the generalization of this method to non-relativistic quantum 
statistics [211 ISl], and, in 1960, Lee and Yang improved the treatment to de- 
scribe, in principle, all the perturbation orders [33j- Finally, in 1969, Dashen, 
Ma and Bernstein extended this method to a relativistic quantum system where 
the interactions are computed through the S matrix [34] . Each one of these 
treatments apply to a different statistical system, but all of them were con- 
structed so as to be valid on a plane static space-time. This prevents their 
straight application to cosmology, where one needs to consider the possibility 
of a curved manifold. This means that the perturbation methods for modeling 
real gases are valid only if the space-section of the universe is plane, and the 
expansion is of the quasi-static type. And these requirements, as said in section 
[21 are fulfilled in the pre- nucleosynthesis universe {kT > 20 MeV) . 

According to the cluster expansion technique, the grand canonical potential 
for a one-component fluid is 

oo oo 

n{z,T) = ^ 6^ = H ^JV e^''('^"«+™) . (27) 

N=l N=l 

The bN are the cluster integrals and encapsulate all the information about the 
interaction processes. The Appendix is a resume on cluster expansions, with 
the differences between classical and quantum systems particularly emphasized. 
The first cluster integrals for the non-relativistic classical system are 

= '^Wi ^^^^^ 

b2 = f{ri,r2)d\^Sr2, (28b) 



^ -9 QX^y J J y [/(^i'^2)/(ri,f3)/(f2,f3) + /(fi,f2)/(ri,f3)-h 

+f{ri,r2)f{r2,r3) + /(ri, f3)/(n^, ^3)] d'ridV2dV3, (28c) 

where A is the thermal wavelength, g is the degeneracy degree and 

/(r;,f,)^e-'^^(^~-^^^)-l, (29) 

are the Mayer functions. Notice that in the classical case the interaction is 
introduced through the intcrparticlc potential V{fi,rj). Though the system is 
nonrelativistic, the rest mass is already included in the 6jv through the factors 
g-/3m^ These are not considered by the traditional texts on non-relativistic Sta- 
tistical Mechanics, but they arc necessary in the derivation of coherent cosmo- 
logical energy density. The pressure and the numerical density arc not affected 
by these factors. 
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Dashen, Ma and Berstein [33] have shown that the general form of the coef- 
ficients b]\r for a RQS is given in terms of the S-matrix as 



where is the cluster integral of the non-interacting quantum theory, A is 
the symmetrization operator, S is the on-shell S-matrix operator [3S], and Cn 
stands for all the A^-particle connected diagrams. We are using the definition 

See the Appendix for further details. 

Once the grand canonical potential p7p is obtained, it is straightforward to 
calculate the pressure ()llap , the numerical density (jllbp and the energy density 
(jllcp as functions of the temperature and the fugacity: 



p{z, kT) ^kTY^ b^z^ , (32a) 

N=l 

oo 

(z,/cT) = ^ iVfeATZ^ , (32b) 



N=l 



pi^,kT) = ikTf±^z-. (32c) 

N=l ^ ' 

These equations are the one-component fluid EOS in parametric form. 

An alternative description is obtained if everything is rewritten in terms of 
the numerical density after inversion of the n(z, T) series. The result is the 
virial expansion, 

oo 

p(n, kT) ^kT^ai {kT)n} , (33a) 
1=1 

oo 

p{n, kT) = (fcT)^ J2 ikT)n^ , (33b) 



1=1 



where a/ and c/ are, respectively, the virial coefficients for the pressure and the 
energy density. They are completely determined by the Bn- For instance, the 
first three terms are: 



U2 



«i = 1- «2 = -|, «3 = |(2^-63j ; (34) 

- J__96]_ _ J_ / nT_ 3\ _ da2 _ 1 ^03 , . 

~ bid (kT) ~ kT\kT 2) ' d{kT) ' " 2 d{kT) ■ ^ ' 
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The question asking for an answer now is: what is the most suitable set of 
EOS for Cosmofogy, Eqs. ((32al[32bl[32cl) or p3a[[33b)) ? The foUowing comments 
wiU help to determine the prcfcrred set. 

The great advantage of writing the EOS in form of series is the possibility of 
approximating them by their first terms. Whenever series are involved, some- 
thing must be said about their convergence. The sets of EOS (|32al I32bl I32cp 
and (|33ai I33bp are equivalent when all the terms of their series are considered; 
but that equivalence ceases to exist when the series are truncated at a given 
order. The convergence of the sets {p{z,T), p{z,T)} and {p{n,T), p{n,T)} de- 
pends critically on the interaction processes at play and, in consequence, on the 
cosmological period under consideration. 

All the relevant particles for the dynamics of the PNU are born through 
pair- production at high temperatures (from one tenth to hundreds of MeV). 
Each of the species will have [U [6j a nearly vanishing (total) chemical potential 
/z ~ and, therefore, a fugacity close to one, z ~ 1. Consequently, z is not a 
good expansion parameter for a series. 

On the other hand, the expansions of p{n, T) and p(n, T) describe with 
great accuracy rarefied non-relativistic gases. Indeed, using the Lennard- Jones 
potential, the correction 

oc 

a2 = -- I (e-^^'--)/^^ - l)rl^dr,2 , = \r2 - ri| , (36) 

9 J 



leads to an EOS that fits the experimental curves p x n for several gases with 
very good precision [3^ . This is a strong physical argument in favor of the set 
{p{'n^T)^ p{n^T)] and indicates that the associated series are well-defined. 

There is another reasoning suggesting the convergence of the series in n: the 
virial coefficients a; and c; are evaluated from a subclass of connected diagrams, 
the irreducible diagrams [lOj . In the classical case, the irreducible diagrams are 
easily detectable, since they are multipy connected^ i.e., each particle connects 
to at least other two. For example, coefHcient 63, given by (|28c|) . is constructed 
from all the four diagrams shown in the right of Figure 2. The coefficients 03 
and C3, however, depend solely on the last diagram of Figure 2. The first can 
be written as 

03 = j I /l2/l3/23rfVi2dVi3 . (37) 

[Both in and ([57)1 we are supposing that the potential depends only on the 
interparticle distance, V{ri,rj) =V{\ri — rj\)-] This restriction on the diagrams 
for the non-relativistic classical system favors the virial expansion against the 
series in the fugacity. It is to be expected that convergence is better achieved 
by the virial expansion. 

In view of these arguments, we choose the set composed by T) and 
p{n, T) as the most convenient for computing interactions in the pre-nucleosynthcsis 
cosmological universe. Nevertheless, even this set presents a validity limit. The 
perturbative methods do not apply to dense systems (with high values of n) or 
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Figure 2: Classical diagrams representing the first three cluster integrals. Each 
line corresponds to a Mayer function and each ball, to a particle. 

to systems submitted to long-range interactions. In such cases the EOS must 
be found by other methods. 

In next section, we will construct a toy-model in order to exhibit some possi- 
ble effects in the evolution of the universe when interaction between the source 
components is taken into account. We shall also argue in favor of the con- 
vergence of the series p{n, T) in this particular example, comparing with the 
behavior of the expansion p(z, T). 

4 Example of the interaction influence 

In our toy-model the primeval interacting fluid is constituted only by photons 
and nucleons coming from pair production. The photons 7 will be treated in 
the usual manner as utra-relativistic ideal bosons, and the nucleons N will be 
considered non-relativistic interacting classical particles. 
The interaction processes taken into account are: 

1. Creation and annihilation of the nucleons in the thermal bath, 7 + 7 ^ 
N -\- N. This reaction generates the mean numerical density of nucleons 
n^v and anti- nucleons n^. In a purely classical context, this process does 
not take place. Here, it serves only as a source of the interacting particles, 
which are treated classically as soon as they come into existence. 

2. Nucleons affecting nucleons through a (very simplified) nuclear potential. 

The electromagnetic interaction cancels out globally since the numerical den- 
sity of nucleons are identical to that of anti- nucleons N (Debye scenery). Weak 
interaction is several order of magnitude less effective than the strong interaction 
and it is consequently neglected. 

For simplicity, we admit that the interactions NN, NN and NN are de- 
scribed by the same nuclear potential (charge independence of the strong inter- 
action), and that the internal degrees of freedom come from spin and isospin. 
Therefore, a hadronic part of the cosmic fluid is composed by particles with mass 
niN = iripf = 938.26 MeV (the proton rest-mass) and degeneracy = (7^ = 4. 
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Figure 3: Potential modeling the nuclear interaction. 



The nuclear interaction shall be modeled by a square-well combined with a 
hard-core potential, as in Figure 3. 

The choice of that nuclear potential is justified by two facts: (i) it presents a 
behavior similar to those exhibited by some successful phenomenological nuclear 
potentials, as those described in Refs. [37l [38j; and (ii) with this simplified 
potential, it is possible to calculate analytically the cluster integrals to the third 
order. 

The parameters of the square-well hard-core potential are obtained from (a) 
the deuteron binding energy, the proton and deutcron mean-squared radius - 
which fix the well's width {b — c) = 1.3 fm and its depth, Vq = 75.6 MeV - 
and (b) from nucleons high-energy scattering data - setting the extension of the 
hard-core c = 0.4 fm (cf. Ref. [39]). 

The pressure and the energy density for the PNU are then written in the 
form: 

p{kT,nN) — pj{kT) + pi\r{kT,nN) , (38a) 
p{kT,nN) = p^{kT) + pN{kT,nN) , (38b) 

where pn and pN are given by (|33aP and (|33bp respectively, and [10] 

p^ TT^kTr 

P7 y ; P-f(kT) = — — . (39) 

Truncating the nucleons EOS in the third order (approximation valid for un 
not too large), explicit expressions forp(T', n^v) and plT^n^) result: 

p{kT, un) ~ ^ + (kT) [un + a2n% + a^nU , (40a) 
45 

p{kT, un) ^ + {rriN + ^kT)nN + {kTf [cjn^ + c^n^ , (40b) 
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Figure 4: Graphics of the pressure p{kT) and energy density p{kT), both mea- 
sured in GeV/fm^, as functions of kT given in GeV. Full lines represent p and 
p of the proposed model (with interaction). For sake of comparison, the dotted 
curves show p and p in the ideal case (without interaction). 



with 

UNikT) ~ 9N^-^^ + 262 + 363 . (41) 

Recall that /ijv — in the pre- nucleosynthesis cosmological period and, in con- 
sequence, Ziv = 1. Coefficients 02 and 03 are found analytically with the help 
of and (P7|) . while &2, ^3 and C2, C3 are derived directly from and ([551) . 
Performing all calculations and inserting njv(fcT) in (|40aP and (|40bp . we finally 
obtain the EOS p(fcT) and p{kT) for the PNU. The result is presented in a 
graphic form - see Figure 4. For kT < 0.3 GeV, the proposed model is quali- 
tatively identical to the ideal case. This is an expected feature since, until 0.3 
GeV, the numerical density n^^kT) is too small to cause any relevant interac- 
tion effect. From this energy value on, the deviation from the ideal case emerges, 
first in the curve for p and then in the plot of p. As energy increases, the two 
functions tend to decrease and, eventually, p and p become negative. This pe- 
culiar characteristic is due to the action of the interaction terms in (|40ap and 
(|40bp . They dominate at high-energies: the attractive part of the nuclear po- 
tential (square-well) make 02 and C2 the most important terms of the truncated 
expansion. 

The effect of the nuclear interaction in the PNU is better seen by observing 
the equation of state in its modern-cosmology familiar form, 

""^^""^ = Tm ■ 

The behavior of w{kT) is shown in Figure 5. For kT < 0.1 GeV the function 
w{kT) 1/3, implying that the EOS is that typical of a radiation-dominated 
universe. The ideal terms reduce slightly this value of w as the energy increases. 
But, for energy values greater than 0.3 GeV, the interaction processes become 
relevant, reducing abruptly the value of w{kT) and deviating its behavior from 
the ideal case. 
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Figure 5: Parametric equation w{kT) as a function of kT (in GeV). The full 
curve represents the proposed model (with interaction). The dotted line corre- 
sponds to the ideal (without interaction) cosmological EOS. 



The dynamical (interaction) terms modify the cosmological EOS and, con- 
sequently, change the form of the primordial expansion. In particular, our 
toy-model produces an accelarated expansion which naturally evolves to a de- 
celerated radiation-like expansion, which is necessary for the nucleosynthesis. 
Nevertheless, we would not dare to say that this simplistic model represents 
realistically the PNU dynamics. This model is apt only to show the importance 
of considering interactions within the cosmic fluid as an important factor in the 
determination of the early universe's evolution. 

Another argument favoring the better convergence of the pressure series 
p{kT, un) (|40ap in comparison to that for the energy-density scries p{kT, zn) 
(|32ap follows as a by-product. It is possible to compare the second interaction 
terms (02 or 62) with the third ones (03 or 63), by defining the functions 



F„(fcT) = 1 - 



0371^ 



F,{kT) = 1 - 



(43) 



The more are Fn{kT) and Fz{kT) close to 1, the more the third interactions 
terms are irrelevant compared to the second terms, and the larger is the pos- 
sibility of convergence of p(kT,nN) and p{kT,zpf). Notice, however, that this 
analysis does not proves convergence: it uses only the first terms of the series. 
Fn{kT) and Fz{kT) only help us to find an argument favoring the good behavior 
of the EOS. Graphics of F„ and Fz as functions of kT are showed in Figure 6. 
The conclusion is that the series in terms of un - namely p{kT,nN) - has a 
better chance to converge than the series in zn- This corroborates the choice of 
section I3.1[ where we have chosen the EOS set given in terms of the numerical 
density as the suitable EOS for the pre-nucleosynthesis cosmological universe. 
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Figure 6: The functions F„{kT) (full line) and F,{kT) (dashed line). kT is 
measured in GeV. 

5 Final remarks 

This work highlights the absence of dynamical terms as direct source of curva- 
ture in the usual cosmological models. The detailed examination of this fact is 
performed in section[51 The analysis is trivial for the equilibrium approach, but, 
even in the non-equilibrium case, no mention is made of the possible relevance 
of interacting components for the universe evolution. Some exceptions are Refs. 

A method for including these dynamical effects in the pre-nucleosynthesis 
universe has been proposed, through the corrections of the EOS based in the 
cluster expansion technique of Statistical Mechanics. The approach is rather 
general and allows the treatment of the cosmic fluid as a classical or quantum 
system (relativistic or not). In principle, the method can be applied to other 
periods of cosmic history, provided the conditions for thermodynamical equilib- 
rium and series convergence are respected^ 

A toy-model has been presented which illustrates the deep consequences 
which interaction between constituents can have for cosmic evolution. Even if 
overmuch simplified, the model points the way toward more realistic approaches, 
based on fundamental physics. More accurate models for the pre-nucleosynthesis 
universe would include other particles than just photons and nucleon-anti- 
nucleons pairs: at least pions, kaons, electrons and neutrinos should be included. 
Besides, the nuclear interaction should be treated in a more complete manner 
than just a square-well hard-core potential. More realistic EOS are found in 
studies [3HI331I31] concerning interacting hadrons. 

We cannot affirm categorically that the inclusion of dynamical terms can 
describe properly the primordial acceleration of the universe (as suggested by 

Such conditions are required only for the interacting part of the cosmic fluid. 
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our toy- model), but the results presented here are indicative that these terms 
should not be simply ignored, as usually done. Maybe the inflationary era 
and the present-day acceleration are just consequences of neglected interaction 
terms. 



A Appendix: Cluster expansions 

The ensemble approach to Statistical Physics is able to include interactions 
perturbatively, via the cluster expansion formalism. The method is developed 
in the grand canonical ensemble, whose partition function, written for a single 
component fluid, is: 



E{z,V,T) = J2Qn{V,T)z 
with z is given by (jlOp and 



N 



N=0 



Qn{V,T) 



W^{ri,...,r]s)<fri...<frN, 



(44) 



(45) 



where tag X indicates the nature of the system under consideration: X = C 
denotes a non-rclativistic classical system and X = Q, a quantum system (QS). 
Wjy is the probability density for a A^-particlc system, 



W^^(fi,...,r-W) 



— (3m 



N 



exp 



N 



i<j=i 



M/^(fi,...,fw) = 7V!(fi,...,fjvMe-'^^"|ri,...,rlv) , 



(46a) 



(46b) 



where Vifiyfj) is the potential between the particles i and j; A is the sym- 
metrization operator and Hj^ the iV-particle Hamiltonian operator. 

The thermodynamical quantities are associated to the grand canonical po- 
tential Q which, just like the partition function, can be expressed as a series in 
terms of the fugacity z. 



-l oo 

f}(z,T) = -lnS(z,y,r)= ^5^z^ 



with 



5Ar 



NIV 



U§in,...,rN)d^ri...(frN 



(47) 



(48) 



The are the Ursell functions and g counts the degeneracy coming from 
internal degrees of freedom|f| The coefficients are the cluster integrals. 

* How to count the degenerate states depends criticaly on the type of system. For instance, 
when deaUng with a quantum system, the counting must respect the symmetry condition of 
the complete wave-function. 
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Using Eqs.lgl]), (gS]), (gT]) and (gS]) and expanding InS in terms of the fu- 
gacity, one writes Uf^ in terms of . Examples: 

U^{n,r2) - M/2^(n,r2) - VFf^(ri)W^i^(r2) , 

W^in,f2,r3) - 3Wi'{n,r2)Wi'ir3) + 2Wi' in)Wi' if2)Wi' (r^) . (49) 

In principle, all the possible states of a statistical system can be decomposed 
in generic diagrams accounting for the several correlation processes (interactions 
or quantum exchange effects). The Qn, determined from the probability density 
functions, are obtained as the sum of all distinct diagrams of the N particles. 
The cluster integrals, determined from the Urscll functions, are obtained as the 
sum of all connected diagrams, 

bjy ^ (connected diagrams of N particles) . (50) 

In a iV-particle connected diagram all the N particles are linked directly or 
indirectly. The links include all kinds of correlation. 

The arrangement of the several distinct diagrams in connected diagrams is 
complicated and it will not be carried out here. The interested reader may 
consult Refs. [T^ for classical systems, and [33] for the quantum case. 

A.l Non-relativistic classical system 

According to (|46ap : 



U^in) = W,^'in) = ^, (51) 



—f3m 



2 



U^irur2) = l^-^j firun), (52) 

where 

/(r--,n,)^e-^^(^^-^^')-l (53) 

is the Mayer function. 

Also from (|46ap . it is possible to show that the non-relativistic classical 
cluster functions can be decomposed in products of Mayer functions [35]. For 
instance: 

(fi,f2,f3) - (e-^^(*'i'*'=)e-^^(^^i'^^^)e-'5^(-~=^^^^) - Se-'^^^*'-^^^) + 2) . 



3 



')( 



A3 

g-/3V(ri,r2)g-/3V(ri,f3) g-/3V(ri ,r2) g-/3y (fs .fa) j_ g-/3V(ri .fj) g-/3y (Fa ,r3) 
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/g-/3m\ 

C/f (ri,r2,r3) = ( j [/(^i: ^2)/(?"i, ?"3)/(?"2, j'a) + /(ri, r2)/(ri, 7%) 

r2)/(r2, rs) + f (n, f^) f (f^ , f^)] . (54) 

Such a decomposition can be obtained for all orders. The Mayer functions will 
always be the fundamental entities in the determination of the [/^(ri, ...,f]\[), 
and lead to a representation of the cluster integrals in terms of connected dia- 
grams — the first three are shown in Figure 2. 

The number of connected diagrams grows very fast as N increases: there is 
only one diagram if = 2, 4 diagrams for = 3, and 38 forA^ ~ 4. 



A. 2 Relativistic quantum system 

Relativistic quantum system is much more complicated than the non-relativistic 
classical system by several reasons: 

• In a relativistic system an interaction cannot be represented by a potential. 
In fact, the determination of an interaction Hamiltonian H' = H — Hq, the 
total Hamiltonian H minus the free Hamiltonian Hq, is a difficult task. 

• There are quantum correlation effects which are mixed to the dynamical 
(interaction) terms. 

• There are discrete bound states. 



Despite these difficulties, there exists a general formalism giving the coeffi- 
cients &Ar in terms of the S-matrix [34] . 

Analogously to (j46bp . we define the operator Un such as 



?7^(fi, ...,nv) = A^! (fi, ...,fjv| Un \ri, ■•■,nv) 



(55) 



Then, using (|48 



5JV 



V 



(fi, ...,fjv| Un In, ...,r7v) (fri...d^rN 



(56) 



or through a Fourier transform. 



Un 



ki, kN) d^ki...(fikN ■ 



The general form of the bN is 



(57) 



bN = yTrUN 



(58) 



Equation (|56p is the realization of this trace in coordinate representation, and 
(|57p is its realization in the momentum representation. 



20 



In order to obtain the &Ar in terms of the S-matrix, it is necessary to separate 
the statistical effects from the dynamical effects, because only the latter influ- 
ence the S-matrix. The traditional method to accomplish this is to construct 

b'^'-^TrU^^\ (59) 

where index (0) denotes a free system (without interaction), t/^'' is related to 
a free system, it is responsible only for statistical contribution. Therefore, it is 
sufficient to subtract ([55)1 from ([55]) to retain the pure dynamical terms. This 
gives 

bN - b^N^ = ^Tr (Um - U^^'>) . (60) 
Thus, the grand canonical potential ([47| is: 

oo 

n-n, = Y,{bN- 5^') , (61) 

N = l 

where flo is the ideal grand canonical potential © • Dashen, Ma and Bernstein 
[34] showed that ([60)) can be written in terms of the S-matrix as: 

Tr{u.-U^^^)-l'-^Tr(A.S-^^s] dE , (62) 

\ / Cjv 

where A is the symmetrization operator, S is the on-shell S-matrix operator 
[35] ■ and cjv indicate that only the connected A^-particle diagrams are to be 
considered. 

Taking ([S^ into (|60|) leads to the general expression (|30p for the cluster 
integrals Bn- That form of bpf includes not only the scattering states, but also 
the (bound state) composed particles (for more details see [34]). 

Let us exhibit the explicit expressions of ([62]) for onc-particlc and two- 
particle systems. 

1. iV = 1: 

Tr (J7i - ^7^"') = Tr (Wi) -Tr (w^"'>) = TrAe'^"^ ~TrAe-P"'°' = . 

(^^) 

The last equality results from Hi = ifj"' (free particle). 

2. iV = 2: 

Tr (U2 - U^^^^ =TrW2~ Tr (w^iW^i) - (Tr#f ^ - Tr (w^f ^VK^'"')) . 
Tr (U2 - ') = Tr {w-^) Tr {w^'^) . 
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-J^Tr^AS^^-S^JdE, (65) 

where H2 is the two-particle Hamihonian operator and S2 is the S-matrix 
operator associated to H2 ■ It is worth to remark that ([55)1 can be rewritten 
in terms of measurable phase-shifts, if we choose the angular momentum 
representation. 
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